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Abstract A search for a charged Higgs boson (H^) 
in tt decays is presented, where one of the top quarks 
decays via t — > followed by — > two jets (cs). 

The other top quark decays to Wb, where the W boson 
then decays into a lepton (e//i) and a neutrino. The 
data were recorded in pp collisions at y/s — 7 TeV by 
the ATLAS detector at the LHC in 2011, and corre- 
spond to an integrated luminosity of 4.7 fb~^. With no 
observation of a signal, 95% confidence level (CL) up- 
per limits are set on the decay branching ratio of top 
quarks to charged Higgs bosons varying between 5% 
and 1% for masses between 90 GeV and 150 GeV, 
assuming B{H+ cs) = 100%. 

PACS 12.60.Fr, 14.65.Ha, 14.80.Da, 14.80.Fd 



1 Introduction 

In the Standard Model (SM), electroweak symmetry 
breaking (EWSB) occurs through a single complex scalar 
doublet field and results in a single physical state, the 
Higgs boson . A particle with characteristics of the 
SM Higgs boson has been discovered by both ATLAS [4] 
and CMS [51]. Beyond the SM, many models have been 
proposed, extending the Higgs sector to explain EWSB. 
The newly discovered boson is compatible with many of 
these models so that discovering its true nature is cru- 
cial to understanding EWSB. Two Higgs-doublet mod- 
els (2HDM) are simple extensions of the SM with 
five observable Higgs bosons, of which two are charged 
{H+ and H') and three are neutral (ft.", 77" and A°). 
The discovery of a charged Higgs boson would be a sig- 
nal for new physics beyond the SM. 

The Minimal Supersymmetric Standard Model 
(MSSM) y is an example of a 2HDM. At tree level, the 
MSSM Higgs sector is determined by two independent 



parameters, which can be taken to be the mass mH+ 
and the ratio of the two Higgs doublet vacuum expec- 
tation values, parametrized by tan/3. In the MSSM, a 
light H+ (defined as < TTit) decays predomi- 

nantly to cs, bbW~^ , and t'^v, with the respective branch- 
ing ratios depending on tan /3 and r7i^+ . For tan /3 < 1, 
cs is an important decay mode with B{H^ — >■ cs) near 
70% [i,[9| for mH± 110 GeV, whereas for tan/3 > 3, 
H+ ^T+v dominates (90%). 

The LEP experimentsplaced lower limits on mjy+ 
in any type-II 2HDM |10[ varying between 75 GeV 
and 91 GeV depending on the assumed decay 

branching ratios for the charged Higgs boson. At the 
Tevatron, searches for charged Higgs bosons have been 
extended to larger values of m^^ + . No evidence for a H'^ 
was found and upper limits were set on the branching 
ratio B{t — >■ H'^b) varying between 10% and 30% for 
a light under the assumption of B{H^ — > cs) = 

100% [il,[i3. 

In this paper, a search for a charged Higgs boson 
produced in tt decays is presented, where one of the top 
quarks decays via t — > H^b with the charged Higgs bo- 
son subsequently decaying to two jets (cs), where again 
a 100% branching fraction is assumed. The other top 
quark decays according to the SM via i — )■ W~b with 
the W boson decaying into a lepton (e//i) and the cor- 
responding neutrino. The signal process therefore has 
the same topology as SM tt decays in the lepton -I- jets 
channel, where one W decays to two jets and the other 
to a lepton and corresponding neutrino, but the invari- 
ant mass of the two jets from the peaks at 'm[j+ . 
The search is performed by comparing the dijet mass 
spectrum in the data with the prediction from SM top- 



^ Charge conjugated processes are also implied throughout 
this paper. 
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quark decays and with the expectation of a top quark 
having a non-zero branching ratio for decay to 



2 Detector description and event samples 

The data used in the analysis were recorded by the 
ATLAS detector in proton-proton (pp) colhsions at a 
centre-of-mass energy of = 7 TeV during the 2011 
data-taking period of the Large Hadron Colhder (LHC) 
[isj . Events were required to pass a high-transverse mo- 
mentum (pt) single-lepton {e/fi) trigger, and to have 
been recorded when aU detector systems critical to muon, 
electron, and jet identification were operational. The 
lepton triggers required in the different data taking pe- 
riods had varying pT thresholds: 20-22 GeV for the 
electron trigger and 18 GeV for the muon trigger. The 
resulting dataset corresponds to an integrated luminos- 
ity of 4.7 fb-i 

The ATLAS detector f2l]| consists of an inner track- 
ing system immersed in a 2 T axial magnetic field pro- 
vided by a thin solenoid; electromagnetic and hadronic 
calorimeters; and a muon spectrometer (MS) embedded 
in a toroidal magnet system. The inner detector track- 
ing system (ID) comprises a silicon pixel detector clos- 
est to the beamline, a silicon microstrip detector, and 
a straw tube transition radiation tracker. The electro- 
magnetic (EM) calorimeters are high-granularity liquid- 
argon sampling calorimeters with lead as the absorber 
material in the barrel and endcap regions, and copper in 
the forward region. The hadronic calorimetry uses two 
different detector technologies. The barrel calorimeter 
(jryl < 1.7)11 consists of scintillator tiles interleaved with 
steel absorber plates. The endcap (1.5 < \ri\ < 3.2) and 
forward (3.1 < \ri\ < 4.9) calorimeters both use liquid 
argon as the active material, and copper and tungsten 
respectively as the absorber. The MS consists of three 
large superconducting toroids each with eight coils, and 
a system of precision tracking and fast trigger cham- 
bers. 

The largest background to the charged Higgs boson 
signal is the SM production and decay of tt pairs. Addi- 
tional background contributions (referred to as non-U 
backgrounds) arise from the production of a single top 
quark, of a or Z boson with additional jets, of QCD 
multi-jets, or of dibosons. 



^ ATLAS uses a right-handed coordinate system with its ori- 
gin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The 2;-axis 
points from the IP to the centre of the LHC ring, and the 
y-sods points upward. Cylindrical coordinates (r, </>) are used 
in the transverse {x, y) plane, 4> being the azimuthal angle 
around the beam pipe. The pseudorapidity is defined in terms 
of the polar angle 6 as r) = — lntan(0/2). 



Top-quark pair and single top-quark events {Wt- 
channel and s-channel) were generated using the mc@n- 
LO 4.01 [H-iil Monte Carlo (MC) generator coupled 
to Herwig 6.520.2 [26| to provide the parton shower- 
ing and hadronization using the AUET2-CT10 [23,|28| 
tune; Jimmy [i^ was used to model the underlying 
event. Single top-quark events in the t-channel were 
generated using AcerMC 3.8 [s^I coupled to Pythia 
with the AUET2-MRST2007LO** 



27 



6.425 [31| 

tune. W/Z+iet and diboson events were generated us- 
ing the leading-order (LO) Alpgen 2.13 |33| genera- 
tor interfaced to Herwig with the AUET2-CTEQ6L1 
27|,|34| tune. The W/Z+jet simulated data include ded- 
icated samples for heavy-fiavour production (bb, cc and 
c). Signal samples of ti H'^bW'b were generated us- 
ing Pythia 6.425 for seven different 7J+ masses from 
90 GeV to 150 GeV. 

The data are affected by the detector response to 
multiple pp interactions occurring in the same or neigh- 
bouring bunch crossings, known as pile-up. Minimum- 
bias interactions generated by Pythia 6.425 [31] were 
overlaid on the simulated signal and background events. 
The events were re-weighted to reproduce the distribu- 
tion of the number of interactions per bunch crossing 

HQ 



IS 



observed in the data. A Geant4 simulation 
used to model the response of the ATLAS detector, and 
the samples are reconstructed and analysed in the same 
way as the data. 



3 Physics objects and event selection 

Jets are reconstructed from topological clusters of calori- 
meter cells [stI using the anti-fct algorithm 38, 3^ with 
a radius parameter R = 0.4. Topological clusters are 
built using an algorithm that suppresses detector noise. 
Jets are corrected back to particle (truth) level using 
calibrations derived from Monte Carlo simulation and 
validated with both test-beam [i^ and collision-data 
studies I41I] . Events are excluded if they contain a high- 
Pt jet that fails quality criteria rejecting detector noise 
and non-collision backgrounds [4^ . Tagging algorithms 
identify jets originating from &-quark decays by select- 
ing jets with tracks from secondary vertices or those 
with a large impact parameter significance. A multi- 
variate algorithm (MVl) 



431, which uses a neural net- 



work to combine the weights from multiple tagging al- 
gorithms, is used to identify jets originating from b- 
quarks. Jets passing the MVl selection are referred to 
as 6-tagged jets. The selection on the discriminating 
variable of the algorithm achieves an average per-jet ef- 
ficiency of 70% to select &-jets in tt events, with a prob- 
ability to incorrectly tag light jets of less than 0.1% 44j . 
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Muons are required to be identified in both the ID 
and MS, and their momentum is obtained through a 
combined fit of all hits in both systems. Muons are also 
required to satisfy isolation criteria to reject those origi- 
nating from heavy-flavour decays and hadrons misiden- 
tified as muons. The sum of the transverse momenta of 
ID tracks within a cone of AR = ^ (^0)^ + {Ar\)^ ~ 
0.3 around the muon, excluding the muon track itself, 
is required to be less than 2.5 GeV. The transverse 
energy measured in the calorimeters within a cone of 
AR = 0.2, excluding the energy associated with the 
muon, is required to be less than 4 GeV. In addition, 
muons are removed if they are found within AR < 0.4 
of a jet that has pt > 25 GeV (isl liel. 

The reconstruction of electron candidates starts from 
a seed cluster in the second layer of the EM calorime- 
ter. The cluster is matched to a track found in the ID 
and a set of selection criteria are applied to reject elec- 
tron candidates originating from jets . Electrons are 
required to be isolated in order to suppress the QCD 
multi-jet background. The calorimeter isolation is per- 
formed using a cone of AR = 0.2 and the track isolation 
uses a cone of radius AR — 0.3. The calorimeter and 
track isolation cut values are chosen to achieve 90% effi- 
ciency with respect to selected electron candidates [i^ . 
As in the case of muons, the electron itself is excluded 
from the sum over the isolation cone. 

Energy deposits in the calorimeter are expressed as 
four- vectors (£',p), where the direction is determined 
from the position of the calorimeter cluster and the 
nominal interaction point {x = y = z = 0). The clus- 
ters are formed assuming E = \p\. The missing trans- 
verse momentum (E™^^^) is given by the negative of the 
vector sum of the calorimeter four-momenta, projected 
into the {x, y) plane. The E^^^^ calculation uses the en- 
ergy scale appropriate for each physics object described 
above. For muons, the momentum measured from the 
combined tracking is used as the energy. The remaining 
calorimeter cells not associated with any physics object 
are included at the electromagnetic energy scale of the 
calorimeter |49|. 

A set of requirements is imposed to select events 
containing tt decays in the lepton-|-jets channel (isj . 
First, events are required to contain a primary vertex 
with at least five associated tracks to suppress non- 
collision backgrounds. Exactly one electron with a large 
transverse energy {E^ > 25 GeV) and \r]\ < 2.5, exclud- 
ing the barrel-endcap transition region 1.37 < I77I < 
1.52, or one muon with large transverse momentum 
{pt > 20 GeV) and |?7| < 2.5 is required. The selected 
lepton must match a lepton trigger object that caused 
the event to be recorded. Jets present in W/Z +iet events 
tend to originate from soft gluon emissions. These back- 



grounds are therefore reduced by requiring at least four 
jets with Pt > 25 GeV and \ri\ < 2.5. At least two 
jets must be identified as originating from a 6-decay us- 
ing the MVl algorithm. To suppress backgrounds from 
QCD multi-jet events, the missing transverse momen- 
tum is required to be i?™'*''' > 20(30) GeV in the muon 
(electron) channel. Further reduction of the multi-jet 
background is achieved by requiring the transverse mas^ 
(mx) of the lepton and E^^^^ to satisfy mx > 30 GeV 
in the electron channel and {E!^^^^ + tot) > 60 GeV in 
the muon channel. These requirements favour the pres- 
ence of a, W boson, decaying to £1/, in the final state. 
The selections are more stringent in the electron chan- 
nel because of the larger multi-jet background. 



4 Kinematic fit 

In the selected events, the two jets originating from 
the decay of the i7+ must be identified in order to re- 
construct the mass. A kinematic fitter 16[ is used to 
identify and reconstruct the mass of dijets from W/H'^ 
candidates, by fully reconstructing the tt system. In the 
kinematic fitter, the lepton, E^^^^ (assumed to be from 
the neutrino), and four jets are assigned to the decay 
particles from the tt system. The longitudinal compo- 
nent of the neutrino momentum is calculated from the 
constraint that the invariant mass of the leptonic W bo- 
son decay products must be the experimental value 
(80.4 GeV) 50]. This leads to two possible solutions 
for this momentum. The fitter also constrains the in- 
variant mass of the two systems (Mv, bjj) to be within 
-Tt = 1.5 GeV of the top-quark mass 172.5 GeV, which 
is consistent with the measured top-quark mass [5l| . 
When assigning jets in the fitter, 5-tagged jets are as- 
sumed to originate from the 6-quarks. The best bbjj 
combination is found by minimizing a for each as- 
signment of jets to quarks and for the choice of solution 
for the longitudinal neutrino momentum, where the five 
highest-pT jets are considered as possible top-quark de- 
cay products. Since the 6-jets are only allowed to be 
assigned to the 6-quarks, and the two untagged jets are 
assigned to quarks from the same charged boson, there 
are two possible jet configurations overall for events 
with four jets, two of which are 6-tagged. For events 
with at least five jets, the two highest-px jets are al- 
ways assumed to be from the top-quark decay products 
(VF/77+ boson or 6-quark) to reduce the combinatorics 
in the fit procedure. The combination with the smallest 



3mT = ^^/2p!^E^^'''^' 



L,rj,±^rj, ^(1 — COS A(j}) where Afji is the azimuthal 
angle between the lepton and the missing transverse momen- 
tum 
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value, Xmin^ is selected as the best assignment. The 
function minimized in the fit is: 



E 

i=^,4jots 



T ~ Pt ) 



E 



SEJ,fit SEJ,mcas\2 
-Pi ) 



(1) 



'SEJ 



E 



k—jjb.biiJ 



tt\AOOO 

I12000 

> 

LU 

10000 
8000 
6000 
4000 



2000 



0, 



ATLAS 

tfi = 7TeV:|Ldt = 4.7 fb"' 
Data 
SM tt 
Non-tt 
^ SM with uncertainty 




10 20 30 40 50 60 70 80 90 



In the first term, the fitted transverse momenta of 
the lepton and the four jets currently under consid- 
eration are allowed to vary around the measured val- 
ues using the corresponding measured resolutions {(Ji). 
The vector sum of the energy of the remaining jets 
(pt > 15 GeV) in the event, labelled SEJ, is accounted 
for in the second term, which also includes all the contri- 
butions to the i?™"**^ by objects not already in the first 
term of the function. Since the SEJ term is dom- 
inated by additional jets, the resolution for this term 
is extrapolated from the jet resolution. The third term 
constrains the hadronic {jjh) and leptonic {Mv) top- 
quark candidates to have a mass close to the top-quark 
mass. 

The x^iij distribution for selected events in the data 
agrees well with the expectation from the simulation 
(see Fig. [T]). Events are required to have Xmin < 10 to 
remove poorly reconstructed tt events. This selection 
has an efficiency of 63% for SM ti events. The fit re- 
sults in a 12 GeV dijet mass resolution, which is a 20- 
30% improvement in the resolution of the dijet mass 
distribution {rrijj), depending on the mass of the boson 
studied, as shown in Fig.[21 After the fit, there is better 
discrimination between the mass peaks of the W boson 
from SM decays of tt and a 110 GeV boson in this 
example. 

Table [1] shows the number of events observed in the 
data and the number of events expected from the SM 
processes after the selection requirements. The SM tt 
entry includes events from both the lepton -I- jets and 
dilepton tt decay modes, where the dilepton events can 
pass the event selection if the events contain additional 
jets and the second lepton is not identified. Good agree- 
ment is observed between the data and the expectation. 
The table also shows the number of signal events ex- 
pected for B{t H'^b) — 10%. The signal prediction 
accounts for acceptance differences due to the different 
kinematics of the t H^h events relative to the SM 
t — Wb events. 



Fig. 1 Comparison of the distribution of Xmin from the 
kinematic fitter for data and the expectation from the back- 
ground estimates for the combined electron and muon chan- 
nels. The MC simulation is normalized to the expectation for 
the SM (H(t H+b) = 0). The uncertainty shown on the 
background estimate is the combination in quadrature of the 
±l(j systematic uncertainties. The final bin also contains the 
overfiow entries. 
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Fig. 2 Comparison of the dijet mass distribution before 
(upper part) and after (lower part) the kinematic fit and the 
< 10 selection criterion. The distribution is shown for MC 
simulations of SM tt decays and the mjj+ = 110 GeV signal 
(tt — ^ H~^bW~b). The curves are normalized to the same 
area. 
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Channel 


Muon 


Electron 


Data 


10107 


5696 


SM tt W+bW-b 
W/Z + jets 

Single top quark + Diboson 
QCD multi-jet 


8700 ±1800 
420 ±120 
370 ±55 
300 ±150 


5000 ±1000 
180 ±47 
210 ±31 
130 ±63 


Total Expected (SM) 


9800 ±1800 


5500 ±1000 




m„+ = 110 GeV | | 


B{t -^H+b) = 10% : 
tt H+bW-b 
tt W+bW-b 


1400 ±280 
7000 ±1400 


800 ±160 
4000 ±800 


Total Expected {B = 10%) 


9500 ±1700 


5300 ±1000 



Table 1 The expected number of events from SM processes, 
integrated over the full range of dijet masses and the observed 
number of events in the data after all the selection require- 
ments. The expected number of events in the case of a signal 
with = 110 GeV and Bit H+b) = 10% is also shown. 

The tt — >• W'^bW~ b numbers include both the lepton + jets 
and dilepton decay channels. The uncertainties are sum of the 
contributions from statistics and systematic uncertainties. 



5 Systematic uncertainties 

The background estimates and the estimate of the sig- 
nal efRciency are subject to a number of systematic un- 
certainties. The QCD muhi-jet background is estimated 
using a data-driven method [52.] that employs a likeli- 
hood fit to the E^'^^^ distribution in the data, using 
a template for the multi-jet background and templates 
from MC simulations for all other processes. The uncer- 
tainty on the QCD multi-jet background is evaluated to 
be 50% by studying the effect of pile-up events on the fit 
results and by performing likelihood fits on the m'Y{W) 
distribution. The dijet mass distribution of multi-jet 
events is obtained from a control region in the data, 
where leptons are required to be semi-isolated, such 
that the transverse momentum of the inner detector 
tracks in a cone of radius AR = 0.3, excluding the lep- 
ton, satisfies 0.1 < p^^^'^-'^ /p'T{e, ^jl) < 0.3. This selec- 
tion criterion is used rather than an anti-isolation cut, 
as a large fraction of the QCD multi-jet events have a 
very small opening angle between the lepton and the as- 
sociated jet, leading to a bias in the reconstructed top- 
quark mass. Leptons in the control region are also re- 
quired to have a large impact parameter with respect to 
the identified primary vertex (0.2 mm < \do\ < 2 mm) 
and an impact parameter significance \do\/(Jdo > 3. 

The rate of VF+jet events is estimated by a data- 
driven method [53| that uses the observed difference in 
the number of and W~ bosons in the data and the 
charge asymmetry (VF+ — W~)/{W'^ + W~), which 
is calculated to good precision by the MC simulation 



of M^+jets events. The heavy flavour fraction of the 
W^+jets MC simulation is calibrated using + 1 jet 
or 14^ + 2 jet events in the data. The uncertainty on 
the W+jets background is 26% (28%) for the electron 
(muon) channel. The shape of the rrijj distribution for 
W^+jets events is obtained from simulation. 

Uncertainties on the modelling of the detector and 
on theory give rise to systematic uncertainties on the 
signal and background rate estimates. The following 
systematic uncertainties are considered: integrated lu- 
minosity (3.9%) [H, |2^, trigger efficiency (3.5%/l% 
for electron/muon), jet energy scale (1-4.6%) [4l|, jet 
energy resolution (up to 16% smearing) 5J|, and 6-jet 
identification efficiency (5-17%). The last three uncer- 
tainties depend on the pt and rj of the jets. Uncer- 
tainties on lepton reconstruction and identification ef- 
ficiency are determined using a tag and pro be method 
in samples of Z boson and J/^/j decays [55|. The mo- 
mentum resolution and scales are determined from fits 



to samples of W boson, Z boson, and J/ip decays [48|, 
[5^. Additional px-dependent uncertainties are placed 
on the 6-jet (up to 2.5%) and c-jet (up to 1.3%) en- 
ergy scales |4lj. Uncertainties on the modelling of the 
tt background are estimated using a second MC gener- 
ator (PowHEG [57-;59j) and comparing the effect of us- 
ing Pythia and Herwig to perform the parton show- 
ering and hadronization. Uncertainties on initial and 
final state radiation (ISR/FSR) are assessed using Ac- 
erMC interfaced to Pythia and examining the effects 
of changing the ISR/FSR parameters in a range con- 
sistent with experimental data [g^l- The predicted SM 
ti cross-section for pp collisions at ^/s = 7 TeV, ob- 
tained from approximate next to next to LO QCD cal- 
culations, is au = 167^Jg pb for a top-quark mass 
of 172.5 GeV [6l[. The uncertainty on the predicted 
value includes the uncertainty in the renormalization 
and factorization scales, parton density functions, and 
the strong coupling constant. An additional uncertainty 
on the tt cross-section (4.5%) is included due to the 
uncertainty on the top-quark mass. The uncertainty 
on the top-quark mass is 0.9 GeV from the combined 
measurement 5l| at the Tevatron. However, this re- 
sult would be biased in the presence of a — ^ cs 
signal in the lepton + jets channel, so a larger uncer- 
tainty of 1.5 GeV is taken, which is consistent with 
the latest top-quark mass measurement in the dilep- 
ton channel from the CMS experiment [H^. Changing 
the top-quark mass leads to altered event kinematics, 
which results in a final uncertainty on the event rate of 
1.9%. The effects of these systematic uncertainties on 
the overall normalization are listed in Table [5] The jet 
energy calibration, 6-jet identification, ti background 
modelling, and ISR/FSR uncertainties also modify the 
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Systematic Source 



Shape dependent 



Jet energy scale 


±9.5% 


6-jet energy scale 


+0.3, -0.6% 


c-jet energy scale 


+0.1, -0.3% 


Jet energy resolution 


±0.9% 


MC generator 


±4.3% 


Parton shower 


±3.1% 


ISR/FSR 


±8.8% 


Shape independent 


b-tagging efficiency (fo-jets) 


±11% 


b-tagging efficiency (c-jets) 


±2.4% 


b mistag rate 


±1.8% 


Lepton identification 


±1.4% 


Lepton reconstruction 


±1.0% 


t-quark mass 


±1.9% 


tt cross-section 


+10, -11% 


Luminosity 


±3.9% 



Table 2 Effect of the systematic uncertainties on the event 
rate of tt background and signal (m^+ = 110 GeV) events be- 
fore any reduction from the likelihood fit, described in Sect.[6l 
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Fig. 3 The dijet mass distribution from data and the ex- 
pectation from the SM {B = 0). The error bars represent the 
statistical uncertainty on the data. The uncertainty shown 
on the background estimate is the combination in quadrature 
of the ±1(7 systematic uncertainties, accounting for the con- 
straint from the profile likelihood fit. The first and last bins 
contain the underfiow and overflow events respectively. 



shape of the dijet mass distribution and are therefore 
determined as a function of rrijj. The systematic un- 
certainties that affect the shape of the mjj distribu- 
tion (top half of Table [5]) are more important than 
the shape-independent uncertainties. The effects of the 
systematic uncertainties are comparable, within 10%, 
between the SM and signal tt samples. The combined 
uncertainty on the single top-quark and diboson back- 
grounds is 15%, which comes mostly from the uncer- 
tainties on the cross-section, jet energy scale, and b- 
tagging. The total uncertainty on the overall normal- 
ization of the non-tt backgrounds is 30%. 



6 Results 

The data are found to be in good agreement with the 
distribution of the dijet mass expected from SM pro- 
cesses (see Fig. [3]). The fractional uncertainty on the 
signal-plus-background model is comparable to the back- 
ground only model. Upper limits on the branching ratio 
B{t -> are extracted as a function of the charged 

Higgs boson mass. The upper limits are calculated as- 
suming the charged Higgs always decays to cs. The fol- 
lowing likelihood function is used to describe the ex- 
pected number of events as a function of the branching 
ratio: 



^(^'") = 11 ri 11 



(2) 



systematic uncertainty. The number of expected signal 
plus background events in each bin, Vi(B,a), is given 

by 



v,{B, a) = 2B{1 - B) auCA"^ Sf^ H pf («J 

i 

-(1 - Bf a^.CA'^Sr n + n pU'^^ 



(3) 



where rii is the number of events observed in bin i 
of the dijet mass distribution and j labels the sources of 



where is the expected number of non-tt back- 
ground events, au is the cross-section for tt produc- 
tion, C is the integrated luminosity, B is the branching 
ratio oi t H^b, and and A^ are the accep- 

tances for signal {ti H+bivb) and SM tt (tt — jjbivb 
and tt — )■ Ivblvb) events respectively. The decay mode 
tt — )■ H^bH^b does not contribute to the expectation 
because this mode does not produce a single isolated 
lepton and hence has a negligible efficiency to pass the 
selection requirements. The S^^ {Sj^) parameter de- 
scribes the shape of the rujj spectrum (normalized to 
one) for 7?+ (W) boson production. It gives the rela- 
tive number of events in bin i according to the normal- 
ized rrijj distribution. The aj variables are nuisance 
parameters representing the systematic uncertainties, 
which are constrained via the Gaussian terms in Eq. [2l 
The Pji functions account for the effect of nuisance 
parameters on the yields and are defined such that 
Pji{aj = ±1) represents the ilcr fractional change in 
the number of entries in bin i of the dijet mass spectrum 
due to systematic uncertainty j. The physics measure- 
ment involves a sufficiently large number of events that 
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Higgs Mass 


Expected limit 


Observed limit 




(stat.® syst.) 


(stat.® syst.) 


90 GeV 


0.080 


0.051 


100 GeV 


0.034 


0.034 


110 GeV 


0.026 


0.025 


120 GeV 


0.021 


0.018 


130 GeV 


0.023 


0.014 


140 GeV 


0.020 


0.013 


150 GeV 


0.015 


0.012 



Table 3 Expected and observed 95% CL limits, including 
systematic uncertainties, on the branching ratio for a top- 
quark to decay to a charged Higgs boson and a 6-quark, as- 
suming that B{H~^ cs) = 100%. The limits shown are 
calculated using the CLg limit-setting procedure. 



CD 0.08 
c 

^0.06 
^ 0.04 
0.02 














Limits at 95% CL: - 


: 


ATLAS 


Expected Limit ^ 




. \!s = 7 TeV 1 Ldt = 4.7 fb 


Expected ± 1a \ 




^ B(H*^c5) = 100% 


Expected + 2o ^ 




—— Observed Limit - 



























90 100 110 120 130 140 150 



m„. [GeV] 

Fig. 4 The extracted 95% CL upper hmits on B{t -> H+b) 
are shown in the range of the charged Higgs mass from 90 
GeV to 150 GeV. The limits shown are calculated using the 
CLg limit-setting procedure. 



this likelihood can constrain the aj parameters beyond 
the precision of the subsidiary measurements. The ef- 
fects of systematic uncertainties are applied coherently 
in signal and background distributions. The subsidiary 
measurements of the aj parameters are taken to be un- 
correlated. The fit uses 17 nuisance parameters in to- 
tal. None of them are shifted by more than one sigma 
compared to the original values obtained in subsidiary 
measurements. Maximal reduction of uncertainty is ob- 
tained for the jet energy scale parameter which is re- 
duced by 50%. 



7 Conclusions 

A search for charged Higgs bosons decaying to cs in 
tt production has been presented. The dijet mass dis- 
tribution is in good agreement with the expectation 
from the SM and limits are set on the branching ra- 
tio B{t H+b), assuming B{H+ cs) = 100%. 
The observed limits range from B — 5% to 1% for 
mH+ = 90 GeV to 150 GeV. These are the best lim- 
its to date on charged Higgs boson production in this 
channel. 



The limits on the branching ratio are extracted us- 
ing the CLs technique at 95% confidence level [63ll64|. 
The consistency of the data with the background model 
can be determined by comparing the value of the test 
statistic (a profile likelihood ratio based on Eq. [2]) in 
the data with the expectation from background-only 
Monte Carlo simulated experiments. The correspond- 
ing probability (p- value) for the background to produce 
the observed mass distribution varies from 67% to 71% 
as a function of TO//+ , indicating that there is no signif- 
icant deviation from the background hypothesis. The 
expected and observed limits, shown in Table [3] and 
Fig. m are calculated using asymptotic formulae (63| . 
The expected limits on B, including the statistical un- 
certainty, vary between 1-3% depending on mjj+ . The 
observed limits, including both statistical and system- 
atic uncertainties, vary between 1-5%. The extracted 
limits are the most stringent to date on the branching 
ratio B{t H+b), assuming B{H+ cs) = 100%. 
These results can be used to set limits for a generic 
scalar charged boson decaying to dijets in top-quark 
decays, as long as the width of the resonance formed is 
less than the experimental dijet resolution of 12 GeV. 
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